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Abstract In this Letter, we argue that the standard analytical derivations of properties of HII regions, such as the speed, shape 
and asymptotic position of ionisation fronts require a more precise treatment. These derivations use the on the spot approxi- 
mation, which in effect ignores the diffuse component of the radiation field. We show that, in fact, HII regions are diffusion 
dominated. This has as a result that the morphology of inhomogeneous HII regions will be drastically different, because shad- 
owing effects have a less profound impact on the apparent shape. Moreover, it will have influence on the propagation speed of 
ionisation fronts. We quantify our claims by analytically deriving the internal radiation structure of HII regions, taking diffusion 
fully into account for several different cosmologically relevant density distributions. 
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1. Introduction 

Ionized hydrogen is produced when the UV radiation emitted 
copiously by hot young stars ionizes the surrounding interstel- 
lar and intergalactic medium. Recently, the resultant HII re- 
gions have gained a lot of interest, because when the first hot 
and massive Pop III stars form, the emitted energetic photons 
blow ionisation bubbles, which after overlapping will eventu- 
ally fill the Universe at the end of the Epoch of Reionisation. 
Therefore, much effort has been put into the physical, analyti- 
' cal and numerical understanding of the mechanisms involved. 
A very relevant property of HII regions, is the structure 
of the diffuse radiation field, as this will influence the dy- 
namics of the expansion, but also the morpholo gy of shadows 
formed within the region. It is well known (Hong & Sung 1989; 
Lopez- Martin etalJf200ll) that diffuse photons dominate near 
the Stromgren radius, but for practical reasons this is often ig- 
nored, under the assumption that it is a small effect. In this 
Letter, we quantitatively show that it is a non-negligible effect, 
especially in the presence of density gradients. 

2. The On-The-Spot Approximation 

Most analytical work done on HII regions heavily depends on 
the on-the-spot approxi mation, part of which has its origin in 
Baker & Menzel ( 1938), and h as b een discussed an d used in 
the standard texts lSpitzetl i ll 9781) and|Osterbrock Jl989l) . In this 
argument, one considers the gas within the ionisation front, or 
within the Stromgren sphere when the static equilibrium so- 
lution has already been reached, and does some bookkeeping 
to equate the number of recombinations within that volume to 



the number of ionising photons. Electron captures directly into 
level n = 1, parametrised by the coefficient a\{T), produce 
photons energetic enough to ionise another atom (the spectrum 
of these diffuse photons can be approximated by a delta func- 
tion just above the Lyman limit). The ionisation balance can be 
drawn up as follows: 



4nr 2 n(r)— = S , + Ana x (T) \ n 2 (r)r z dr 
at J 
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a(T) J n 



2 (r)r 2 dr 
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in which r is the position of the front and S * is the number of 
ionising photons (i.e. with frequencies v > vo above the Lyman 
limit treshold) emitted by the source per second. The second 
term on the right-hand side is the total number of diffuse pho- 
tons, given the HI density distribution n(r). The last term on 
the right-hand side accounts for the total number of recombi- 
nations to every level, parametrised by the recombination coef- 
ficient oa(T). The integrations are over the whole HII region. If 
the Stromgren radius has already been reached, the three con- 
tributions on the rhs of Eqn.Q cancel, but up until that point 
the ionisation front has a finite speed dr/dt. 

Unlike the source photons, which are directed outwards, 
the diffuse photons may cross the nebula in any direction from 
their point of creation. Thus, one has to do a full treatment of 
radiative transfer, which complicates matters drastically. This 
is where the on-the-spot approximation enters the picture. For 
typical HII regions, the physical parameters are such that the 
opti cal depth for photons near the Lyman limit is ~ 30 (cf. 
IOsterbrocMfl9 89). Thus, one argues that the diffuse photons 
produced by recombinations to the ground level will be re- 
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absorbed very close to where they were produced (on the spot). 
This has as a result that effectively the recombinations directly 
to the ground level do not count, because they are exactly bal- 
anced by the photons they produce. Henceforth, one can ignore 
these recombinations, and the diffuse photons created by them, 
and use the effective C ase B (using terminology introduced by 
baker & Menzell938l) recombination coefficient 



a B (T) = a A (T) - ai (T), 
with Eqn.Q reducing to 

n dr 
4nr 2 n(r) — = S t + 47ra B (T) 
at 



(r)r 2 dr. 



(2) 



(3) 



The Stromgren sphere radius can be easily obtained from 
Eqn.Q by considering the equilibrium solution with zero ve- 
locity of the ionisation front drjdt = 0. 

We will argue that it is not permitted to use the on the spot 
approximation for anything more than just calculating the ra- 
dius of the Stromgren sphere. First and foremost, the on the 
spot approximation is based on the argument that the mean free 
path for Lyman limit photons is very small, and that thus all dif- 
fuse photons are locally re-absorbed. But this argument is not 
correct, because the mean free path for the source photons (as- 
suming these are also near the Lyman limit) is just as small. 
Thus, the diffuse and the source photons are on equal footing 
energy-wise, and no distinction between them can be made, ex- 
cluding their directionality. Thus, at each location r within the 
HII region, not only are the diffuse photons used to compen- 
sate for the a\ recombinations, but also a weighted fraction of 
the source photons, by which a fraction of the diffuse photons 
can actually escape. In effect, a fraction of the mono-directional 
source photons is converted into diffuse radiation. 

One can therefore see that the number of directional source 
photons decreases, when moving outwards from the source, un- 
til at a certain point the diffuse radiation starts to dominate. In 
the next section, we will analytically derive the diffuse versus 
mono-directional structure of the radiation within a Stromgren 
sphere for several astrophysically important density distribu- 
tions. 



3. The diffuse structure within Stromgren spheres 

In this section, we will derive analytically the diffuse and 
source radiation structure within a spherically symmetric HII 
region which is already in equilibrium, i.e. the Stromgren ra- 
dius has already been reached. In subsequent sections, we will 
use these results to draw conclusions for more general cases. 



3.1. Homogeneous HI matter distribution 

In the following, we assume a homogeneous hydrogen medium 
density n m (r) = n Hh a temperature T = 10 4 K for the ionised 
plasma, and, following the considerations that led up to the on 
the spot approximation, we assume that the mean free path of 
the photons is very small, until locally the gas is fully ionised. 
We do the calculations for a Stromgren sphere in d-dimensional 



space, because we can use the general result in the following 
subsections. 

Because the Stromgren radius has already been reached, 
we know that there is an equilibrium between the number of 
recombinations per timestep and the total number of ionising 
photons emitted by the source and the gas, so we may drop the 
time-dependence. The number of recombinations to every level 
and the number of diffuse photons produced within a spherical 
shell at distance r is 



N mc (r) 
N m (r) 



f(</)ar A nV -1 </r 
Z(d)a l n 2 r d - 1 dr, 



(4) 



in which £ = In ^ 2 /F(d/2) (= 4tt in 3D) is just a geometrical 
factor and T(x) is the Gamma function. The number of photons 
emitted by the source plus the number of diffuse photons within 
the whole region must compensate all the recombinations, from 
which we obtain 



= %(d)a B 



n 2 r d - 
Jo 



'dr 



-—a B n R s , 

d & 



(5) 



in which R$ is the Stromgren radius. 

If we define 7 s (r) and /<j(r) as the number of source and dif- 
fuse photons per unit volume, respectively, left at radius r < R&, 
we can write down the following coupled system of differential 
equations on the domain < r < R& 



ffdteAnV- 1 /.^) (7,(r) + 7 d (r))" 



2„d-\ 



dr 

= -f(<0a A BV- 1 / d (r)(/ s (r) + / d (r))-' +£( £ /)a 1 «V 

The first term on the rhs of both equations is the total number 
of recombinations being compensated for via a weighing term, 
that takes care of the fact that, if there are more source than dif- 
fuse photons present, the recombinations will mainly be com- 
pensated by the source photons, and vice versa. This system 
of equations is closed via the initial conditions /d(0) = and 
7 S (0) = S ». We should note that we have used a symmetry con- 
dition here, in the sense that the diffuse radiation sent into the 
opposite direction (away from the front) is exactly balanced by 
the diffuse radiation on the other side of the source. 

If we add the two equations, the weighing terms add up 
to unity, and using the initial conditions, we easily obtain an 
equation for the total radiation field 



(/.+/ d )(r)= ^a B n 2 R d s 



(7) 



We can plug this into the weighing terms in Eqs.(|6j, and upon 
integrating we obtain 



h(r) = (/, + 7 d )(r) - 7 d (r) = ^a B n 2 R d s 



'"'5 



■(8) 



After choosing units such that ^■a B n 2 Rj - 1, and using 
the ratio a A /a B = 4.18/2.60 at 10 4 K (e.g. lOsterbrocW 989). 
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Figure 1. Graphs of the functions (I s +Id)(r) (solid), I s (dashed- 
dotted), Id (dashed), and Id/h (dotted) for a homogeneous 
(top), an r _1 (centre), and an r~ 2 (bottom) medium distribution. 



we can draw a figure that plots (/ s + Id)(r), h(r), Id(r) and the 
ratio Id(r)/I s (r). This has been done for the normal d = 3 space 
in Fig.^ top. One can see that at a certain radius r eq the diffuse 
radiation component takes over. We can derive r eq exactly by 
solving I s (r) = Id(r), and obtain 

r eq =R s (l-2- a »' a >f ld . (9) 
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Filling in the values results for 3D in r eq /7?s = 87.6%. 
Thus, we can immediately conclude that for a homogeneous 
gas distribution the outer 12% of the HII region is dominated 
by diffuse radiation, and that very near the front only diffuse 
photons are present. This means that, in terms of volume, 33% 
of the HII region is diffusion dominated. This effect is even 
greater for higher temperatures, for which the ratio a B /a\ is 
smaller. 

3.2. An r~ { HI gas distribution 

The effect described in the previous section is even more drastic 
for a HI density distribution in the form n(r) = n c (r/r c )~ l . This 
density profile has been used to de scribe the central regions of 
a dark matter halo jNavarro et alJl997l) . 

Following the same arguments as in the previous 
Subsection, we get for the 3D case 

( M-ecO) = 4na A n 2 .r 2 dr 

\ N diB (r) = 4na in y c dr (10) 
[S. = 4na B ny c R s , 

which are similar to Eqs.@ and Q for d = 1, but with a differ- 
ent constant factor. Thus, we can use Eqs.Q, (|SJl and and 
obtain 

' (/, + Id)(r) = 47ia B n 2 r 2 R s (l - 

</,(/•) = 47ia B ny c R s (l-f s )^ (11) 

The graphs of these functions are plotted in the centre part 
of Fig.H Filling in the value for a B ja\ at T — 10 4 K, gives 
r eq /Rs - 67.2%, from which we can conclude that for an 
isotropic r~' distributed HII region the outer 33% is diffusion 
dominated, which is equivalent to roughly 70% of the total vol- 
ume. 

3.3. An r~ 2 HI gas distribution 

The effect is strongest for a HI density distribution in the form 
n(r) = n c (r/r c )~ 2 , which has, for example, been used to de- 
scribe the density profile of halos just interior to the virial ra- 
dius and of collapsing isothermal spheres. In the derivation 
for this density distribution, we have to be more careful, be- 
cause the equivalents of Eqs.Q and (|5} will contain singulari- 
ties at r — 0, which can be resolved by choosing an inner edge 
< ro < 7?s- Thus, we obtain 

N iec (r) = 4na A n^r~ 2 dr 
< N m (r) = \na x ny c r- 2 dr (12) 
S. = 4na B n 2 y(r l - R s l ) . 

Using similar steps as in Subsection 3.1, we obtain 

' (7, + I d )(r) = 4na B ny c (r~ l - R & 1 ) 

- h(r) = 47ia B n 2 c r A c (r Q l - R^f (r" 1 -R?)* (13) 

r eq = R s (l+2- a *'<*(R s /ro-l))~ 1 . 

The graphs of these functions are plotted in the bottom part of 
Fig-Dior ro < r < R$. Taking a typical value of ro/^s = 0.05, 
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r/R_s 

Figure 2. Graphs of the functions (Id/I s )(r) for the varying val- 
ues (from top to bottom) c = 1.0, 0.9, 0.8, 0.7, 0.6, 0.5 in a r~ 2 
density distribution. Overplotted is the line Id/h - 1 which 
determines the radius above which the diffusion dominates. 

we find that r eq /7?s = 13.8%, which means that, in terms of 
volume, 99.7% of the HII region is diffusion dominated. 

3.4. Hard photons 

The analysis in the previous subsections was done assuming 
that the source spectrum is peaked just above the Lyman limit. 
It is much more realistic to assume that a substantial part of the 
luminosity will be in the higher frequency domain. These pho- 
tons have a longer mean free path (a(v) °c v~ 3 ) than the diffuse 
photons, from which we expect that their contribution will in 
effect move r eq towards Rs. To accurately model this, we would 
have to do a full radiative transfer treatment, but we can give 
an analytical estimate of this effect. Hereto, we do not take the 
source spectrum to be a delta function just above the Lyman 
limit, but at a certain higher frequency. As a result, we relax 
the condition that the medium is optically very thick to source 
photons, by which not 7 s (r) photons contribute to the locally 
available radiation field, but only a factor < c < 1 thereof. 
Nothing changes in the equations except that the weighing fac- 
tors will now be c/ s /(c/ s + Id) and Id/(cl s + Id). We solved the 
equations for the r~ 2 distribution for varying values of c and 
plotted the results in Fig. |2] One can see that the position of 
r eq does change, but that the diffuse radiation field effect is still 
very noticeable. 

4. Concluding 

What does all of this analysis entail? First and foremost, it 
shows that a large part of HII regions, in some cases even al- 
most the whole volume, is diffusion dominated. This has a pro- 
found effect on the formation of shadows behind clumps within 
the HII region. A more quantitative analysis on this matter has 
to be done in forthcoming work, but one can immediately see 



that, even in the case of the only moderately diffuse homoge- 
nous medium distribution (subsection 3.1), if a dense clump is 
put within the radius r eq , the resulting shadow will be squeezed 
by the diffuse field in the outer region. This ef fect w ill be much 
stronger than in the analysis in Canto etalJ dl998h . in which 
only the locally produced diffuse photons can move into the 
shadow region, whereas in our analysis also the overall dif- 
fuse radiation field will have a profound influence. This effect 
will have an impact on (cosmological) radiative transfer codes 
which use ray tracing methods to delineate sharp shadows, and 
the resultant 21cm signatures of the HI regions in between the 
overlapping (re)ionisation bubbles, which may be observed in 
the near future. 

Another effect is that the diffuse radiation field will influ- 
ence the dynamics of the HII region. From the considerations 
in the previous sections for the equilibrium solution, one can 
see that the only radiation component pushing against the front 
is the diffuse one, which is not mono-directional, but isotropic. 
Thus, the net flux ionising photons in the direction of the front 
is less. The symmetry condition, as used in subsection 3.1, is 
only valid in the equilibrium (time-independent) solution, but 
in the time-dependent case the speed of photons diffusing from 
one side of the HII region to the other is finite, and, as a re- 
sult, this effect may slow down the ionisation front, by which 
the standard a nalytical solutions to Eqn.Q, as found in e.g. 
Spitzejl lll978l) . will have to be modified to take into account 
that this diffusion speed is considerably smaller than the speed 
of light. If it turns out that this effect can not be neglected, it 
will have an influence on the analysis of how fast the first stars 
and/or quasars could have reionised the Universe. 

Concluding, the usual on-the-spot approximation in which 
the diffuse radiation component is neglected can only be used to 
calculate the Stromgren sphere within an isotropic HI medium 
distribution. We have shown that the diffuse component can 
not be neglected, which will result in a drastic change of the 
morphology of anistropic distributions, and which may change 
the speed of the ionisation fronts. 
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